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METHOD AND SYSTEM FOR CONCURRENT ESTIMATION OF FREQUENCY 
OFFSET AND MODULATION INDEX 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This patent application claims the benefit of priority from and incorporates by 
reference the entire disclosure of U.S. Provisional Patent Application No. 60/392,114, filed 
on June 27, 2002. 

5 

BACKGROUND OF THE INVENTION 

Technical Field of the Invention 

The present invention relates generally to the field of radio receivers that utilize 
continuous phase modulation (CPM) and, more particularly, to a method of and system for 
10 estimating a modulation index and a carrier frequency offset of a CPM signal 

Description of Related Art 

Wireless technologies, such as terrestrial and satellite mobile communications and/or 

BLUETOOTH systems, may use continuous-phase-modulated (CPM) signals to transmit 
15 data. Binary CPM or M-ary CPM may be employed for the wireless transmission of data 

packets. When data is transmitted using CPM, the modulation index may need to be known 

in some receiver architectures. 

To improve performance, receiver architectures could be employed that require 

knowledge of the value of the modulation index of the transmitted signal. Due to the use of 
20 independent frequency generating circuits in the transmitting and receiving devices, a carrier 

frequency offset is typically generated. In order to achieve optimal performance, the carrier 
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frequency offset should be compensated for as much as possible. Therefore, there is a need 
for a method of and system for estimating the modulation index and the carrier frequency 
offset of a CPM signal. 

5 SUMMARY OF THE INVENTION 

These and other drawbacks are overcome by embodiments of the present invention, 
which provide a method of and system for concurrent estimation of a modulation index and 
frequency offset of a CPM signal An estimator for estimating a modulation index and 
frequency offset of a received continuous-phase-modulated (CPM) signal includes at lea,st 
10 two filters for filtering the received CPM signal, a calculator for calculating an a value and a 
(3 value, and a processor for receiving a signal output by each of the at least two filters, the a 
value, and the (3 value. The processor is adapted to calculate estimates of the modulation, 
index and frequency offset from the signals received by the processor and the received a 
value and (3 value. 

1 5 A method of estimating a modulation index and frequency offset of a received 

continuous-phase-modulated (CPM) signal includes filtering the received CPM signal via at 
least two filters, calculating an a value and a (3 value, receiving a signal output by each of the 
at least two, filters, the a value, and the P value, and calculating estimates of the modulation 
index and frequency offset from the received signals and the received a value and P value. 

20 An estimator for estimating a modulation index and frequency offset of a received 

continuousrphase-modulated (CPM) signal includes a noise whitener for whitening noise of 
the received CPM signal, at least two filters for filtering the noise-whitened CPM signal, an 
initializer for processing a training sequence, and a processor for receiving a signal output by 
each of the at least two filters and the processed training sequence. The processor is adapted 
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to calculate estimates of the modulation index and frequency offset from the signals received 
by the processor and the processed training sequence. 

A method of estimating a modulation index and frequency offset of a received 
continuous-phase-modulated (CPM) signal includes whitening noise of the received CPM 
5 signal, filtering the noise- whitened CPM signal via at least two filters, processing a training 
sequence, receiving a signal output by each of the at least two filters and the processed 
training sequence, and calculating estimates of the modulation index and frequency offset 
from the received signals and the processed training sequence. 

An estimator for estimating a modulation index and frequency offset of a received 
10 continuous-phase-modulated (CPM) signal includes at least two filters for filtering the 

received CPM signal, a noise whitener for whitening noise of a signal output by the at least 
two filters, an initializer for processing a training sequence, a processor for receiving signals 
output by the noise whitener and the processed training sequence. The processor is adapted 
to calculate an estimate of the modulation index and the frequency offset from the received 
1 5 signals and the processed training sequence. 

An estimator for estimating a modulation index and frequency offset of a received 
continuous-phase- modulated (CPM) signal includes a receiver for receiving the CPM signal 
and a processor for estimating the modulation index arid frequency offset according to the 

following equation v = (B T C~ l B)~ l B T C~ l <fi . v represents a vector including elements 
20 representing scaled versions of estimates of the modulation index and the frequency offset. C 
represents a noise covariance matrix, B represents a data model matrix, and (f> is an 
observation vector that represents a phase of the CPM signal. 
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Further advantages and specific details of the present invention will become apparent 
hereinafter from the detailed description given below in conjunction with the following 
drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

A more complete understanding of exemplary embodiments of the present invention 
can be achieved by reference to the following Detailed Description of Exemplary 
Embodiments of the Invention when taken in conjunction with the accompanying Drawings, 
wherein: 

FIGURE 1 is a block diagram that schematically illustrates an estimator in accordance 
with principles of the present invention; 

FIGURE 2 is a block diagram of the estimator of FIGURE 1, including a bias removal 
component in accordance with principles of the present invention; 

FIGURE 3 is a block diagram that schematically illustrates another estimator in 
accordance with principles of the present invention; 

FIGURE 3 A is a block diagram that schematically illustrates an estimator in 
accordance with principles of the present invention; 

FIGURE 3B is a block diagram that schematically illustrates an estimator in 
accordance with principles of the present invention; and 

FIGURE 4 is a block diagram that schematically illustrates a noise-whitening 
estimator in accordance with principles of the present invention. 
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DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS OF THE 
INVENTION 

In accordance with principles of the present invention, several approaches can be 
5 taken to estimate the modulation index and carrier frequency offset of a CPM signal. The 
estimators associated with aspects of the present invention may be divided into two distinct 
categories, namely estimators that assume white noise and estimators that assume colored 
noise. Another classification within the two distinct categories may be made based upon 
whether or not Inter-Symbol Interference (ISI) is assumed to be present in the input signal. 
10 When it is assumed that ISI is present, two further subclasses may be introduced based on 
whether or not the value of a parameter e is known or unknown. The estimators described 
hereinafter are based upon the least-squares approach and can be represented by the 
following equation. 

v = (B 7 C: ] By i B T C l 0 (1) 

15 

In equation (1) above, v is a vector that includes elements representing scaled versions of the 
estimates of the modulation index and the frequency offset. The matrix C represents a noise 
covariance matrix and the matrix B represents the data model. The last three terms in 
equation (1) are a filter operation on an observation vector <f> which is the phase input to the 
20 estimator. 

A first approach uses an estimator based upon a simple data model, which does not 
take into account the Inter-Symbol Interference (ISI). The first approach, as shown in detail 
in FIGURE 1, assumes white noise and no ISI. By assuming that bit timing and frame 
synchronization are known, ah output (i.e., an element of the observation vector) (|> k of a 1- 
25 bit differential demodulator at an optimal sampling moment can be modeled as: 
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(2) 



where bk is a transmitted bit, h is a modulation index, f is an actual frequency offset, T sym is a 
symbol period, and nk is a distortion term that includes noise and ISI. By omitting the noise 
terms nk , equation (2) can be rewritten in matrix notation as follows: 



(f>- Bv in which </> = 



ft>N 



, v = 



and 5 = 



ft; 1 
K 1 



K 1 



(3) 



If (i.e., N transmitted bits) are known, equation (3) can be solved by 

multiplying the observation vector § by the pseudo-inverse of B. Assuming white noise is 

present, the covariance matrix C is equal to the identity matrix. Therefore, equation (1) may 
be rewritten as: 

v = (B T By l B T 0 (4) 
wherein the superscript T denotes the transpose of the matrix B and the superscript -1 denotes 
the inverse of the resulting matrix shown in parentheses. 

Given the above structure of the data matrix B, the following equations prove true: 



B T B 



with S - ^b k , a - 



N S 
S N 

N 



and (B T By ] = 



a 6 
B a 



(5) 



and 0 = 



-S 



k=l 



N 2 ~S 2 ^ N 2 ~S 2 
The value of S directly relates to a weight Wof the known bit sequence as follows: 

S-2W-N (6) 
By implementing a least-squares method, the values of x and y of the vector v can be 
estimated. The estimates of the modulation index h and the frequency offset f may be 
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directly derived from the estimates of x and y respectively. The estimates of the values of x 
and y of the vector v can be obtained by applying two linear operations on the observation 
vector <f) and a post-processing step that depends on the weight of the training sequence. As 
such, equation (4) remains true for v as follows: 



v = 



a 






j 







with q - B T (/) = 



N 

2>* 



(7) 



wherein q x and q 2 are the elements of the vector q , with q x representing the output of the 

first filter (correlation operation) and q 2 representing the output of the second low-pass 
comb filter. FIGURE 1 illustrates a block diagram of an estimator 100 in accordance with 
principles of the present invention. The estimator 100 of FIGURE 1 is based on the data 
10 matrix shown in equation (3). The estimator 100 implements the operations of equation (7) 
and, as mentioned above, assumes white noise and no ISI. A received signal ^ (the signal 

received, mixed down to base-band, and differentially demodulated) is passed through a first 
finite-impulse-response (FIR) filter 102 to yield qj. The coefficients for the correlation filter 
102 are +1 or -1. The received signal <f> k is also passed through a second FIR filter 104 to 
15 yield q2. 

As described above, the training sequence, which is a data sequence known at both 
the receiver and the transmitter, can be used to derive a and /?. The calculated a is output to a 
first multiplier 106 A and a fourth multiplier 106D. The derived /? is output to a second 
multiplier 106B and a third multiplier 106C. qj, which is output by the first FIR filter 102, is 
20 multiplied, at the first multiplier 106 A, with the derived a. q) is also multiplied, at the second 
multiplier 106B, with the derived /?. 
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q 2 , which is output by the second FIR filter 104, is multiplied with the derived /? at the 
third multiplier 106C. q 2 , which is output by the second FIR filter 104, is also multiplied with 
the derived a at the fourth multiplier 1 06D. 

The result output by the first multiplier 106 A and the result output by the third 
5 multiplier 106C are added at a first adder 108 A. The result output by the second multiplier 
106B and the result output by the fourth multiplier 106D are added at a second adder 108B. 
The result output by the first adder 108 A is x from equation (2). From equation (2), x can be 
scaled to yield an estimate of the modulation index h. As shown in equation (2), by 
multiplying x with the modulation index h is produced. Therefore, at multiplier 1 1 OA, x 
10 is multiplied with 1/rc, thereby yielding an estimate of the modulation index h. As is also 
evident from equation (2), the output of the second adder 108B, y can be multiplied by 
l/(27cT sym ) at multiplier 1 1 0B to produce an estimate of the frequency offset f. 

Because of the data model used, the simple estimator 100 might yield biased results. 
For example, a non-zero mean noise term or correlation between the noise and the desired 
1 5 signal might cause the simple estimator 100 to produce unsatisfactory results. The bias in the 
estimate of the modulation index typically depends on one or more of the frequency offset, 
the modulation index, and the value of a signal-to-noise ratio E b /N 0 Most typically, no 
significant bias is present for the estimate of the frequency offset. 

In the event that E b jN Q is difficult to estimate, the bias in the modulation-index 

20 estimate can be compensated for at a particular value of E b /N Q . For example, the value of 
E b /N Q at which the receiver operates at a bit-error rate (BER) of 10" 3 can be selected. In 
addition, because the bias in the modulation-index estimate depends on the modulation index 
itself, the bias can be compensated for at a typical modulation index value, such as, for 
example, 0.32 
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A post-processing step in accordance with principles of the present invention takes 
into account the fact that the bias in the estimate of the modulation index h depends 
approximately quadratically on the estimated frequency offset f in order to compensate for 
the bias in the modulation index estimate. For a bias-compensated modulation index hcomp, 
5 the following quadratic equation holds: 

Kom P =h + c 0 +c 2 y 2 (8) 

The coefficients c 0 and c 2 are chosen via a curve-fitting process in order to minimize the bias. 
The bias in the modulation-index estimate and the frequency-offset estimate can be derived 
by simulation. From the simulation results, adequate bias-reduction processes could be 

10 derived via curve fitting.. 

FIGURE 2 illustrates the estimator of FIGURE 1 with additional post-processing to 
remove bias from the estimate of the modulation index h. The estimator 200 of FIGURE 2 is 
similar, to the estimator 100 of FIGURE 1, except for the implementation of additional 
components used to introduce the coefficients c 0 and c 2 . As noted above, the coefficients c 0 

1 5 and c 2 are used to remove bias from the estimate of the modulation index h. 

As indicated in equation (8), y ? which is output by the second adder 108B, is squared 
by a squaring block 202: An output of the squaring block 202 is multiplied with the value of 
c 2 at a multiplier 204. An output of the multiplier 204 is added, at an adder 206, to c 0 and to 
the estimate of the modulation index h As noted above, the estimate of the modulation index 

20 h is output by the multiplier 11 OA. The adder 206 outputs the bias-compensated modulation 
index h camp . 

The simple estimators 100 and 200 represent relatively computationally-efficient 
implementations; however, simplification of the data model implemented by the estimators 
100 and 200 might not always produce optimal results. Therefore, an estimator based on a 
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more-complex data model than that used in the estimators 100 and 200 can be utilized in 
another embodiment of the present invention. 

FIGURE 3 is a block diagram that schematically illustrates a more-complex estimator 
300 in accordance with principles of the present invention. The more-complex estimator 
assumes white noise in a manner similar to that of estimators 100 and 200. However, the 
more-complex estimator 300 assumes that I SI is present in the signal and further that the 
parameter e (see data model from equation (9) shown below) is unknown. 

In the estimator 300, instead of processing the input signal as modeled by equation (3) 
and implemented in the estimators 100 and 200, the model on which the estimator 300 is 
based is a more complicated model, namely equation (12) shown below. To limit the 
complexity of the estimator 300, a relatively simple ISI model has been assumed. However, 
other ISI models can be used without departing from principles of the present invention. 

In the estimator 300, a linear relationship between a differential output phase 0k of a 
transmitter and the input phase (|) k of the estimator is assumed. The linear relationship of the 
differential output phase 6 k and the input phase (J) k can be modeled as follows: 

^=^_ 1+ (l-2^+^ +1 (9) 

Equation (9) shows the relationship between the input phase <J) k and the output phase 
6k and allows for ISI by the parameter e . If no ISI is present, then the parameter s has a 
value of zero. 

The differential output phase 0* can be represented by the following equation: 

\ =<Pk-<Pk-x = b M 0°) 

wherein <pk is the phase of a transmitted symbol a^. Combining equations (9) and (10) and 
adding the impact of the frequency offset f and the distortion term n k yields: 

<f> k = b k hn - 2b k + b M )dm + InfT^ + n k (11) 
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with c k ^'(b^ x -2l> k +b k+l ) 
Equation (11) shows that the ISI exhibits a relationship with the foregoing bit and the 
following bit However, in more severe ISI cases, the ISI may exhibit a relationship with the 
previous two bits and the following two bits. If the distortion term nk is omitted, equation 
(11) can be written in matrix notation as follows: 









~b 2 




f 










X 






c 3 


1 




(f>-Bv in which (f> = 




, v = 


y 


and B = 


K 




1 


(12) 








z 


















_ b N-\ 


<7H 


1 





When the ISI is not neglected, as mentioned above with respect to equation (11), the 
foregoing bit and the following bit are required, and therefore the index of equation (12) 
10 begins with b 2 and ends with b N ^ . 

If the N=2 transmitted bits {b2...b N -i} are known, equation (12) can be solved by 
multiplying the observation vector (f> with the pseudo-inverse of B y such that 
v = (B T By l B T <p , as shown in equation (4) above. 

Referring again to FIGURE 3, the estimator 300, which is described mathematically 
15 in equations (4), (1 1), and (12), requires, in addition to the filtering and correlation shown in 
the estimator 200, another filter, graphically represented as a middle filter 306. The filtering 
and correlation of filters 302 and 304 operate in a manner similar to filters 102 and 104 of 
FIGURE 2. The middle filter 306 has N-2 coefficients Ck. For the coefficients Ck, the 
following holds: Ck e {0, ±2, ±4}. 
20 Variables x and z are manipulated by multipliers 31 OA and 31 0B in a manner similar 

to that shown for x and y in FIGURE 2 to yield the modulation index h and the frequency 
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offset f In cases where the ISI is dominated by the transmitter characteristics and the 
receiver filter chain, the value of the parameter s from equation (9) may be assumed to be 
known. The parameter e is deduced given the, overall filter chain in the transmitter and 
receive parts of the transceiver. Therefore, the estimator 300 may be simplified by assuming, 
in addition to white noise and ISI, that the parameter 8 is known. Due to this fact, the 
estimator 300 can be simplified and the number of filters utilized reduced as shown in 
FIGURES 3 A and 3B. 

In a first simplified implementation, shown in FIGURE 3 A, the estimator 300A 
includes a correlator filter that is slightly more complex because the filter coefficients are no 
longer +1 or -1 as in the simple estimators 100, 200, and 300. For the first simplified 
implementation of the estimator 300 A, equation (10) is substituted into equation (9), and 
equation (1 1) may be rewritten as follows: 



<|>k ^(sb^ +(l-2e)b k +8b k+1 )h7T + 27ifT sym +n k 
A =d k x + y + n k 

with d k = (8b k _j + (1 - 2e)b k -f eb k+1 ), x = hn and y = 27rfT 5 



(13) 



sym 



15 In the implementation shown in equation (13), the value of the parameter 8 is assumed to be 
known. By omitting the distortion term n k , equation (13) can be rewritten in matrix form as 
follows: 



<f>- Bv in which <j> 



V = 



and B 



(H) 



The implementation of the estimator 300A derived from equation (14) requires two filters: 1) 
a low-pass filter (304) similar to that in the estimators 100 and 200; and 2) a correlation filter 
(302) that is matched to the channel (i.e., a matched correlator). Therefore, the N-2 filter 
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coefficients dk are no longer +1 or -1 but take values of the set {±1, ±(l-2e), ±(l-4e)}. As 
such, the correlation filter is more complex than that of the estimators 100, 200, and 300. 
The modulation index h and the frequency offset f are calculated in a manner similar to that 
of FIGURE 3. The vector-matrix multiplier 308 outputs variables x and y, which are in turn 
5 manipulated by multipliers 3 1 OA and 31 0B to form the estimates of the modulation index h 
and frequency offset f 

The first implementation of the estimator 3 00 A requires the matched correlator The 
matched correlator has increased computational complexity; therefore, a second 
implementation of the estimator 300 with reduced computation complexity, is described 
10 below. 

Referring now to FIGURE 3B, the second implementation of the estimator of 
FIGURE 3 is illustrated. The second implementation is not as complex as the first 
implementation; however, the second implementation is more complex than the estimators 
100 or 200. In a manner comparable to that of the estimator 200, the second implementation 
15 includes a post-processing step that needs to be executed only once. 

; The description of the estimator 300 of FIGURE 3 from equation (12) and (3) 
produces the following equation: 

v - {B T B)~\B T <f> = (B T By l w with w = 

p, q, and r are the outputs of the three filters 302, 306, and 304 as shown in the estimator 300. 
20 If the value of 8 is known, then q, which is the output from the middle filter 306, is not 
required. As shown in FIGURE 3B, the middle filter 306 has been eliminated. With P y 

representing the element of the matrix B T B on row i and column j and because B T B is 
symmetric, the following equation results: 
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1 23 
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£ \2 1 22 r 23 

P P P 

.^13 1 23 ^33 , 



hn 

£fl7T 



(16) 



In accordance with equation (16), the following equations prove true: 



hn = 



and 



Cn^33 ^13 ) (^13^23 ^12^33)^ 

9WT -P n P + P»rHP n P n -P u P 2 i)&x 
- p P _ P 2 — 

1 11 33 M3 



(17) 



(18) 



5 Equations (17) and (18) can be manipulated via a processor 320 to yield the estimate of the 
modulation index h and the frequency offset estimate f. The second implementation of the 
estimator 3 00B provides significant advantages over the estimators 100 and 200 when ISI is 
present. 

A third class of estimators includes noise whitening to further improve the 
10 performance of the modulation index h and the frequency offset f estimators. A specific 
configuration of a noise-whitening estimator is given in FIGURE 4 which is a block diagram 
of a noise whitening estimator. Due to differential demodulation preceding the estimation, 
the distortion term no longer exhibits typical white noise characteristics. Once the 
co variance of the matrix of the noise is known, the estimation process can be improved. 
15 The one-sided autocorrelation function R™ of the noise process n k can be 

approximated by the following: 

K n = [1 -0-5] (19) 
Taking into account the noise covariance matrix C, which may be directly deduced from the 
one-sided autocorrelation function R™, equation (15) may be replaced by equation (1). The 
20 noise whitening is performed by multiplication of the inverse C" 1 of the noise covariance 
matrix C. The matrix multiplication of B T C~ X with the observation vector <f> can be 
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implemented in two ways. In a first option, which is used in the estimator of FIGURE 4, this 
operation is implemented by applying n (n being equal to the number of columns of matrix B) 
filters in parallel (multiplying with B T ). In a second option, the operation above is instead 
performed by two subsequent filter operations, where the first filter operates on the 
5 observation vector <j> to whiten the noise present in that vector, i.e. by multiplication of C" 1 . 
Then the output of this filter is fed to the n filters in parallel (n being equal to the number of 
columns in matrix B) y i.e. multiplication by B T . Both options are functionally the same. In 
the second option, the whitening of the noise is performed explicitly, while in the first option, 
the noise whitening is implicitly performed. The estimator described by equation (1) does not 

10 restrict the values of the filter coefficients to +1 or -1, thereby increasing complexity over 
both the estimators 100, 200, and 300. The estimator described by equation (1) is an 
improved noise-whitening estimator that outperforms the estimators 100, 200, and 300 at the 
cost of increased complexity. 

To reduce complexity of the estimator described by equation (1), a first option is to 

15 quantize the inverse of the noise covariance matrix C. Although the complexity might be 
reduced, the quantization introduces a performance loss in the estimator. A second option is 
to adapt the structure of the inverse of the noise covariance matrix C. A finite-impulse- 
response (FIR) filter may be utilized to whiten the noise. Due to the differential demodulated 
estimator input signal the noise has a high-pass characteristic. Approximated whitening can 

20 be achieved by passing the signal through a low-pass filter. One attractive solution would be 
to use a K-tap comb filter. K may be chosen such that a good balance is obtained between 
performance loss (compared to ideal whitening) and complexity reduction. A third option 
would be to implement the approximated whitening operation by means of a low-pass infinite 
impulse response (IIR) filter. 
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Any of the three options discussed above can be utilized to reduce the complexity of 
the noise whitening filter. The best option should be chosen according to applicable system 
characteristics. 

The principles of noise whitening may be applied to any of the above-mentioned 
5 estimators. Depending on the data model, matrix B, that has been assumed, the estimators 
100, 200, 300, 300A, and 300B may include colored-noise compensation. For example, by 
substituting the data model from equation (3) into equation (1), an estimator similar to that of 
estimator 100 is obtained, except that the estimator is now colored noise compensated. 
Substituting the data model from equation (11) into equation (1) yields a noise- whitened 
10 estimator similar to that of estimator 300. All of the estimators 100, 200, 300, 3 00 A, and 
300B may be altered, by changing the noise model used, to yield colored^noise-compensated 
estimators. 

Referring again to FIGURE 4, the incoming signal <(> is passed to each of a first FIR 
filter 404 and a second FIR filter 406 in order to be low-pass filtered. The first and second 

15 FIR filters 404, 406 operate in a manner similar to those of FIGURE 3B, except that 
A - B T C" X . The FIR filters 404, 406 implicitly whiten the noise based on the values of the 
matrix A. An output p of the first FIR filter 404 and an output r of the second FIR filter 406 
are similar to the p and r values of the estimator 300B, except for the addition of the noise 
whitening. The outputs p and r are utilized in further calculations made in a post processor 

20 420. 

An initialization unit 422 of the estimator 400 receives the training sequence 
{b x ...b n }. Depending on the assumed data model, parameter s may also be required (see 
equations (13) and (14)). Once the matrix B, which represents the data model, is calculated, 
then matrix B T C~ ] Bis calculated. The matrix B T C~ ] B , along with the outputs p and r is 
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passed to the post processor 420 as an initialized training sequence and estimates the 
modulation index h and the frequency offset fin accordance with the above equations. 

While exemplary embodiments of the present invention have been described, it should 
be recognized that the invention can be varied in many ways without departing therefrom. 
5 Although the present invention has been described primarily as being used in, for example, an 
ad hoc wireless system operating according to BLUETOOTH, embodiments of the invention 
can also be used in other systems that utilize CPM. Because the invention can be varied in 
numerous ways, it should be understood that the invention should be limited only insofar as is 
required by the scope of the following claims. 

10 
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